sistema de E.D.0O. 
lineares de primeira ordem 


META: 

Descrever a teoria preliminar sobre sistemas de equações diferen- 
ciais lineares de primeira ordem. 

OBJETIVOS: 

Ao fim da aula os alunos deverão ser capazes de: 

Reconhecer um sistema de equações lineares de primeira ordem. 
Transformar uma equação diferencial de ordem superior em um 
sistema de equações de primeira ordem. 

Identificar uma solução geral de um sistema de equações lineares 
de primeira ordem. 

Verificar as hipóteses do teorema de existência e unicidade. 
PRÉ-REQUISITOS 

Resolução de sistemas lineares. Além dos conhecimentos das aulas 


1,2,6,7e8. 
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Sistema de E.D.O. lineares de primeira ordem 
13.1 Introdução 


Caros alunos, nesta aula começaremos a estudar nosso último con- 
teúdo, abordaremos os fundamentos teóricos acerca dos sistemas 
de equações lineares de primeira ordem. Veremos que os resulta- 
dos abordados aqui nada mais são do que uma generalização dos 


resultados vistos na aula 6. 


13.2 Sistema de equações lineares de primeira 


ordem: Fundamentos teóricos 
Definição 13.1. Uma E.D.O. da forma 
Y' = A(x)Y + B(s), 


onde Y = (y,:-:,yn) ER" ex é a variável independente real é 
dita sistema de equações diferenciais ordinárias linear de 


primeira ordem. Quando B(x) = 0 o sistema acima se reduz a 
eo e 


Este sistema reduzido é conhecido por sistema de equações 
diferenciais ordinárias lineares homogêneo de primeira or- 


dem. Caso contrário, o sistema é dito não homogêneo. 


Considere o sistema de equações diferenciais ordinárias linear não 


homogêneo na sua forma vetorial 
Y' = A(x)Y + B(a). (13.69) 


Além da forma vetorial, podemos escrever o sistema (13.69) na sua 
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forma matricial, dada por 


A au(x) aa(x) - an(x) y bi(x) 

wo an(x) az(x) e: an() y2 bo(a) 
= + 

Yh ani(x) analx) eo ann(z) Yn bala) 


ou na forma 


mn = au(z)y + a(x)yo + + arn(x)yn + bi(x) 
yw — an(x)y + avs(x)yo +: + aon(2)yn + do(a) 
U = Ani(x)yn + An2(1)y2 pers ea AnnlE)Ya + bn(a). 


Exemplo 13.1. Uma E.D.O. linear de ordem n 
an(o)yl) + ---ao(x)y = g(2), 


onde an(x) £ O pode ser vista como um sistema de equações de 
primeira ordem. De fato, primeiramente, escrevamos a E.D.O. 


acima na forma 


1 
y O) [-an-i(x)y ao(x)y + g(a)] 
Chame ny =yyp=Y,y= 4", Un =y"-D. Assim, 
w = ED) 
w = 3 
Yn-1 = Yn 
do = glana(oytD =... ao(a)y + g(a)) 
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Matricialmente, o sistema acima tem a forma 


0 j 0 is 0 ; 
y w 
k 0 0 1 ça 0 
Ya . . . . ya 
* |= ; : : + 
: 0 
0 0 1 
Um ao(z) ala) an-a(2)  anca(a) Yn g(x) 


an(z) ante) an(z) an(z) 


E, vetorialmente, é dado por 


0 1 0 0 

0 0 1 0 
onde Y = (yr ---yn)T,A(x) = : 

0 ai 0 1 

ao(z) a(o) o an-a(m) anil) 


an(z) an(x) : an(x) an(z) 


e B(x) = (0 ---Og(x))”. 


Observação 13.1. 1) A notação usada acima B = ()”, se refere 
a matriz transposta. 

2) Observe que uma E.D.O. pode ser vista como um sistema uni- 
dimensional, ou seja, a matriz A(x) é uma matriz de ordem 1 x 1. 
3)Uma das grandes utilidades de se trabalhar com sistemas está 
na redução da ordem das equações a serem estudadas. De fato, no 
exemplo acima vimos que uma E.D.O. de ordem n pode ser estu- 
dada como um sistema de equações de ordem um, o que facilita 
muito sua análise e aumenta o número de equações que podem ser 


analisadas analiticamente. 


Mãos a obra 


Escreva a seguinte equação diferencial 


y" + In(2)y! + 3xy — 22,2 E (0,00) 


(212) 


Equações Diferenciais Ordinárias AU LA 


13 


em forma de sistema de equações diferenciais de primeira ordem. 


Definição 13.2. Um vetor solução do sistema (13.69) num inter- 
valo 1 é qualquer matriz coluna 
gula) 
V = 
ato) 
cujos elementos são funções diferenciáveis que satisfazem o sistema 


(13.69). 


Exemplo 13.2. Verifique que 


são soluções do sistema 
re yY. 
5 3 
Pela definição, para provarmos que Y e Y são soluções do sistema 
dado, basta-nos verificarmos se ele satisfaz o sistema. 


Observe que 


si —2e"2* = 18eº7 
e; e dad 
2e2% 30e8” 
e 
de-2% 13 gem 
2e22 5 3 “e | 
18e8z 13 3e8% 
308? 5 3 5ebz 
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O problema 
Y'= A(a)Y + Bla), Y(x9) = Yo, 


onderelYEeRPeY=(y ml éuma matriz constante, é 
dito P.V.I.. 

Analogamente a equações diferenciais de ordem superior, sistema 
de equações de primeira ordem também possui um teorema de 


existência e unicidade, o qual enunciamos a seguir. 


Teorema 13.1. (Existência e Unicidade) Suponhamos que os 
elementos das matrizes A(x) e B(x) sejam funções contínuas num 
intervalo comum I que contenha xo. Então, existe uma única 


solução para o P.V.I. 
Y'= A(o)Y + B(2), Y(x0) = Yo, 
no intervalo T. 


À seguir enunciaremos uma série de resultados análogos aos resul- 
tados vistos na aula 6 para equações lineares de ordem superior. 
Esses resultados são fundamentais para nosso embasamento teórico 


sobre resolução de sistemas de equações lineares de primeira ordem. 


Teorema 13.2. (Princípio da superposição )Seja Y1,:--,Y; 


um conjunto de vetores solução do sistema de ordem n homogêneo 
Prue AR (13.70) 
no intervalo 1. Então, a combinação linear 
Y=-qnh+na 


também é solução de (13.70) em 1. 
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Teorema 13.3. Sejam Yi(x) = (yu(2) ya(x) « yua(a)” , Yo(x) — 1 3 
Cyra(a) goal) + gnala))O, eos Yalx) = (yin(a) gan(a) eo ymnl))”, 

n vetores solução do sistema homogêneo (13.70) no intervalo T. 

Então, o conjunto de vetores solução será L.Il. em 1 se, e somente 


se, o Wronskiano 


yn(a) yia(s) ah Yin(ã) 
aC A DAR A 
Yni(a) Ynalx) a al) 


para todo x ET. 


Observação 13.2. Na verdade para que constatemos que o con- 


junto Y,:--,Y é L.l., basta que verifiquemos que 
WA, Mm) £O 


para algum 19 € T. 


Mãos a obra 
Decida se o conjunto de vetores Yi(x) — (cosz e)T Ya) — 


(cosz sen2x)” é LI. ou não. 


Definição 13.3. (Conjunto fundamental)Um conjunto funda- 
mental para o sistema homogêneo de ordem n (13.70) é um con- 


junto de n vetores solução de (13.70) L.I. 


Teorema 13.4. (Solução geral )Seja Y1,--:,Yn um conjunto fun- 
damental do sistema homogêneo (13.70) em IT. Então, a solução 


geral do sistema homogêneo (13.70) no intervalo T é 


Ylv)=ch ++, 
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onde c1,:-- ,Cn São constantes quaisquer e, uma solução geral para 


o sistema não homogêneo (13.69) em I é dada por 
Y()=Y()+Ylv) =ahn ++ + Ya), 
onde Yp(x) é uma solução particular do sistema não homogêneo 


(13.69) em T. 


13.3 Sistemas de equações lineares de primeira 
ordem homogêneo com coeficientes cons- 


tantes. 


Lembremos que no exemplo anterior vimos que uma solução do 


sistema 
VS Ay, 
13 
onde 4 = é dada por 
5 3 
=9 
Vs o = à EE = Ke 
—e2% —1 i 
1 
onde K = e À = —2. Podemos nos perguntar se toda 
—1l 


solução será dessa forma. a resposta é sim. 
Suponha que 
Y=Ke”, 


onde K é uma matriz de ordem n x 1 com coeficientes constantes, 


seja solução do sistema homogêneo 


Y'=AY, (13.71) 
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onde Y ER” e A é uma matriz com coeficientes constantes de 
ordem n x n. (De agora por diante usaremos a notação A para 


representar uma matriz com coeficientes constantes.) Assim, 
Y=A4Y o RD =AKeN es (AK-AR =0, 


(Por abuso de notação usaremos o símbolo 0 para denotarmos uma 
matriz nula.) Como eX? £ O para todo x, podemos dividir a ex- 


A 


pressão acima por e” e obter 


(AK —-)K) =06(AK -»dK)=06(A-AdK =, 


onde Id é a matriz identidade. 
Dessa maneira, para acharmos uma solução não trivial do sistema 
(13.71) da forma 

V=KeM, 


onde K não é a matriz nula, temos que resolver o sistema 
(A-AIDK =0. (13.72) 


Da álgebra linear, sabemos que para o sistema homogêneo (13.72) 


ter solução não trivial é necessário que 
det(A — Ald) = 0, (13.73) 


uma vez que se det(A — Ald) O segue que A — AId possui matriz 
inversa, B, por exemplo e assim, B(A —- AId)K = IdK = 0, 
portanto K seria a matriz nula. 

A equação (13.73) é chamada equação característica da matriz 
A, suas soluções, À, são os autovalores de 4. Uma solução K £ 0 


de (13.72) correspondente ao autovalor À é dita autovetor. 
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Na resolução de (13.73) examinaremos três casos: autovalores reais 
distintos, autovalores reais repetidos e autovalores complexos. Come- 
cemos pelo primeiro caso. 

Autovalores reais distintos. 

Suponha qua a matriz A, de ordem n x n, tenha n autovalores 
reais distintos Ay,---, An, então um conjunto L.I. de autovetores 


K1,:--, K, poderá sempre ser obtido e assim 
YN=KeMt... Y= Kent 


será uma conjunto fundamental para o sistema (13.71) para x € 
(— 00,00). Dessa maneira, a solução geral do sistema homogêneo 


(13.71) é dada por 
Val) = Ge pude Roe ME (—00,00). 


Exemplo 13.3. Resolva o sistema homogêneo Y' = AY, onde 


0-0. 4 
A=| 0 -3 0 
7000 


Pelo que foi exposto acima, precisamos achar os autovalores da 


matriz 4. Assim, temos que resolver a equação 
det(A — Ald) = 0, 
que”, neste caso, é dada por 


—À 0 1 
O -3-A 0 =0S0+DA+-FD(A-2)=0. 
2 2 1-A 
Portanto, os autovalores da matriz À são reais e distintos, a saber, 


A =-3,A, = —l e Às = 2. Então, usando um pouco de álgebra 
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linear, é possível mostrar que existe três correspondentes autove- 1 3 
tores L.l., K1, Ko e K3. Vamos encontrá-los? 
Para À = —3, a equação (A — A Id)K4 = 0 assume a forma 

3 01 KI 0 

0:40: ) ro [=| 0 |; 

22 4 K3 O) 


onde Ki = (K1 K92 83)”. Assim, para encontrarmos o autovetor Ki, 


basta-nos resolver o sistema 


351 + ra =0 


281 + 2K9 + 4k3 — 0. 


Com poucos cálculos achamos a solução 81 = K1,42 = 51,3 = 
—3k1, onde K; é uma constante qualquer. Portanto, escolhendo 
Kk1 = 1, segue que o autovetor, K1, correspondente ao autovalor 
A — —3 é dado por Ki = (15 — 3)”. De maneira análoga cal- 
culamos os outros dois autovetores Ko e Ks3, obtendo os valores 
K - (20 —-297,K = (102). Portanto, a solução geral do 


sistema dado é 


onde c4, c2,c3 são constantes arbitrárias. 


Autovalores reais repetidos 
Suponha agora que a matriz À, de ordem n x n, possui um au- 


tovalor, digamos, A com multiplicidade m, m < n. Neste caso, 


podemos ter duas situações: 


(2 19), 
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1) Para algumas matrizes A é possível encontrar m autovetores 
LI. K1,---, Km correspondentes ao autovalor Ay. Nesse caso, a 


solução geral do sistema contém a combinação 


Re E eos eres dog e 


2)Para outras matrizas A podemos encontrar apenas um autovetor 
K4 correspondente ao autovalor A de multiplicidade m. Neste 


caso, as m soluções L.I. são da forma 


X = Ane” 
Y% = KoyxeMt + KogeM 

m-—l A m—2 A a 
Im = Kmitaye O + mota O + + Kmmedo, 


(13.74) 
onde K;; são vetores de ordem n x 1. Veremos como encontrar tais 


vetores no exemplo a seguir. 


Exemplo 13.4. Resolva o sistema Y' = AY, onde 


10 0 
A=|[22 1 
RR 0 
Resolvendo a equação característica det(A — AId) = 0, obtemos 
(A— 1) = 0. Assim, a matriz A possui um autovalor A = 1 com 
multiplicidade 3. Calculemos o autovetor associado. 
Pelo que apresentamos anteriormente, para encontrar o autovetor, 


K1, devemos resolver a equação 


(A-— MId)K,=0, 
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onde Ki, = (1 &2 &3)”. Dessa maneira, temos 


00 0 Ka O) 
21 +11 no |=| 0 |; 
01 —1 K3 O) 


ou escrevendo de outra forma 


2681 +69 — 63=0 


R9— R3 — 0. 
Resolvendo o sistema acima, obtemos Rj = 0,62 = k2,h3 = ko. 


Ássim, 


Note que obtemos apenas um autovetor associado a A1, a saber, 
qualquer um da forma (0 «2 82)”, onde «2 é uma constante ar- 
bitrária. Observe quaisquer dois autovetores nesse formato serão 
linearmente dependentes. Queremos aqui fazer uma pausa para a 


seguinte observação. 


Observação 13.3. Suponha que temos um autovalor A, com mul- 
tiplicidade 2 e que temos achado, ao resolver o sistema, apenas, por 
exemplo, a seguinte relação K1 — K2 + K3 = O. Com essa relação, a 
solução do sistema seria 81 = &2 — k3,k9 = k9 € kg = k3. Dessa 


maneira, teríamos 


R2 — K3 1 —1 
K = [5 =| 1 |no+ 0 k3 
R3 O) 1 


Portanto, o espaço dos autovetores associados a Ay é bidimensional 


e podemos encontrar dois autovetores L.I. associados ao autovalor 
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AM. Fazendo «2 = 1, «3 — 0 obtemos Ki = (11 0)” e escolhendo 
Ko = 0, K3 — 1 obtemos K» = (-1 01)”. Note que Ki, K» são 


vetores L.I. 


Bom, voltando ao nosso exemplo, escolhendo «9 = 1, temos K, = 


(0 1 1)”. Assim, uma solução para o sistema dado é 
Yi(x) = Ke = Ke” = (0,1,1)7e”. 
Por (13.74), uma segunda solução, Y>, L.I. com Y, tem a forma 
Y = Kijxe” + KoeM?, 


Assim, substituindo essa expressão de Y5 no sistema Y' = AY e 


simplificando, obtemos 


(AK1 — AM Ki)xeM? + (AKo — MKo — Ke” =0. 


Ajx Ay 


Como as funções ve!” e e” são L.l. essa igualdade é válida 


quando 
(A-XMId)K, = 0 


(A-MIdK, = Ki. 
Assim, para encontrarmos a expressão de Y5, basta-nos resolver 
a segunda equação do sistema acima e encontrar o valor de Ko. 
Seja Ko — (k1 &9 K3)”, então a segunda equação do sistema acima, 


pode ser escrita da seguinte maneira, 


2rm1i +ro— n9=1 


ko — ka = 1. 


A solução desse sistema é dada por »1 = 0,89 = k9,483 = 69 — 1. 


Assim, escolhendo «2 = 5, temos K» = (0 5 4)”. Dessa maneira, 
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Y> é dada por 
Y=(011)xel+(054)7e”. 
Agora vamos a terceira solução, Y3, a qual por (13.74) é dada por 
22 
YJ= 6 LAR + Kage? + Ke, 


Assim, derivando essa expressão de Y3 e substituindo no sistema 


dado, encontramos o sistema 


(A-MIJK = 0 
(A-MId)K, = K 
(A-MId)Ks = Ko. 


Ássim, para acharmos a expressão de Y3, basta-nos resolver a ter- 
ceira equação do sistema acima, a qual pode ser escrita da seguinte 


maneira 


2rni +Hkho—kK3=5 


Ko— kg =4. 


A solução do sistema acima é dada por &j = 1/2,K2 = kK92,k3 = 
ko — 4. Assim, escolhendo «2 = 0, temos K3 = (5 O — 4)”, Por- 
tanto, 


2 
1 
= oe +(05 4970 +(50 —47e. 


Logo, a solução geral do sistema dado é 


Y(x) = ch + coYa + ca)s. 
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Observação 13.4. Caso o autovalor À; tenha multiplicidade maior 


do que 3, a fórmula dada acima 


(A-MID)K = 0 
(A-MID)K, = K 


(A-MIDKs = K 
continua seguindo a mesma construção, ou seja, a próxima equação 


seria 


(A MId)Ks4 = K3 
e assim por diante. 


Autovalores complexos 


Se algum dos autovalores da matriz A de ordem nx for complexo, 
digamos, A = a + i8, sabemos que o seu conjugado,A» = a — 
18, também será autovalor de A. Nesse caso, teremos o seguinte 


resultado 


Teorema 13.5. Seja 4 a matriz de coeficientes constantes reais 
de (13.71) e seja K1 o autovetor correspondente ao autovalor com- 


plexo, A = a +15, da matriz A. Então, 
KjeN£ e Ke, 
onde K1 é a matriz conjugada de K1, são soluções de (13.71). 


Assim, pelo resultado acima, uma vez que a matriz A tenha au- 
tovalor complexo, esse autovalor contribuirá com uma combinação 
de funções complexas para a solução geral do sistema. 

Analogamente ao que fizemos com equações lineares de primeira 


ordem quando se tinha raízes complexas, trabalharemos a parte 
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complexa da solução geral até obtermos uma versão que envolva 1 3 
apenas funções reais. Desse modo, observe que, usando a fórmula 


de Euler (e'” = cos0 + isen 0), temos 


Ke = KyeLeilz = K,e““[cos Bx + isen Ba] 


Ke = Ke Ê — Ke““[cos 8x — isen Bx]. 
Pelo Teorema 13.2 segue que 


1 eras ES 
Y = g(Krede + Ke?) 


Yo e? 5l-Kiede + K,e?] 


são também soluções de (13.71), as quais são L.I. em (—00,00). U- 


sando a fórmula de Euler nas expressões de Y e Y5 acima, obtemos 


que 

Y, = e“”[B; cos Bx — Basen B4] 
e 

Y, = e“”[Bo cos Bx + Bisen Ba], 
onde 


1 Est ) Ear 
Bj = ala + Ki), Bo= al Ha + Kil. 


Portanto, a solução geral de (13.71), neste caso, contém a combi- 


nação real de funções 


crie" [By cos Ba — Basen Bx] + coeº” [Bo cos Ba + Bysen Ba]. 
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13.4 Conclusão 


Na aula de hoje, concluímos que os sistemas de equações li- 
neares nos ajuda a resolver equações diferenciais lineares de ordem 
superior. Nessa aula foi apresentado o método de resolução de 
sistemas de equações lineares homogêneas com coeficientes cons- 


tantes. 
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RESUMO 


Na aula de hoje aprendemos como resolver um sistema de equações 


lineares homogêneo com coeficientes constantes. Vimos que tudo 
se resume a analisarmos as raízes de uma determinada equação, 
a qual é chamada equação característica da matriz associada ao 
sistema homogêneo. Nesse caso, podemos ter três opções: raízes 


reais distinats, raízes reais repetidas ou raízes complexas. 


PRÓXIMA AULA 


Em nossa próxima aula continuaremos a resolver sistemas de equações 
lineares de primeira ordem só que, desta vez, trataremos dos sis- 


temas não homogêneos. 


ATIVIDADES 


Atividade. 13.1. Resolva o sistema 
TA = ay + by, yo =by+ay,abeR. 


Atividade. 13.2. Resolva 


y 1 7/10 un 
ww |= 4 —-2 —&4 y> 
y3 —7 —5 —6 ya 


Atividade. 13.3. Resolva o P.V.I. 


e — 30 +2y+4+4z 
cu = dx tda se(0)=2y(D)=0,20)=1, 
E = 4+2y+3z 


a 
/ ul 
(228) 
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Atividade. 13.4. Resolva 
E = 61x-—y 


U- 542. 
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